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We use nonrelativistic supersymmetry, mainly Darboux transformations of the general (one- 



Qj , parameter) type, for the quantum oscillator thermodynamie actions. Interesting Darboux gen- 

U-i ■ 

' eralizations of the fundamental Planck and pure vacuum cases are discussed in some detail with 
relevant plots. It is shown that the one-parameter Darboux-transformed Thermodynamics refers 

^^ ' to superpositions of boson and fermion excitations of positive and negative absolute temperature, 

> : respectively Recent results of Arnaud, Chusseau, and Philippe regarding a single mode osciUa- 

^^ , tor Carnot cycle are extended in the same Darboux perspective. We also conjecture a Darboux 

>^ generalization of the fluctuation-dissipation theorem. 

!>■ ' 

O 



I. INTRODUCTION 



-y \ Recently, Arnaud, Chusseau, and Philippe (hereafter ACP) |I| studied the work done per 

c^ I Carnot cycle by a single mode oscillator of the ideal LC type operating between two baths 

^' of different temperatures. They used the old (1906) prescription of Einstein of discrete huj 

exchanges with the baths. Within this approach, they confirmed the Carnot theory of the 
efficiency of cyclic engines for this case. As known, the energy of the oscillator is given by 
rN ' Planck's distribution 



Up{u,f3) = — + ■ 1 

2 exp[ijhuj) — 1 

One can define an action function 

fp{x) = -^ = - + — _— = -coth — , 2 

oj 2 exp(ax) — 1 2 2 

where x = (3uj. This action plays the role of a generalized force in the process of frequency 
change. In order to evaluate the efficiency, ACP introduced a 'two-variable entropy' 

s(x, y) = xfp{y) - J fp{z)dz = — coth(y ) - In sinh(— ) + C , (3) 
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with the property s{x,x) = s{x), the latter being the usual entropy. Employing standard 
thermodynamic formalism, ACP calculated the efficiency of both reversible and nonreversible 
Carnot cycles. 

ACP made the interesting remark that the oscillator action /p(x) fulfills as particular 
solution the following Riccati equation 

As a matter of fact, the constant vacuum action fy = ^ is also a solution of the same 
Riccati equation, whereas the pure thermal action /t = (nx)-! ^^ ^ particular solution of 
the following equation 

| + «/ + f^0. (5) 

which is a Bernoulli equation (or a particular type of Riccati equation). Here we would like 
to study the thermodynamic consequences of using methods belonging to supersymmetric 
quantum mechanics focusing on one-parameter Darboux transformations of thermody- 
namic actions. 



II. GENERAL SOLUTION OF THE SUPERSYMMETRIC PARTNER RICCATI 

EQUATION 



In Witten's supersymmetric quantum mechanics 0, which is a simple application of 
Darboux transformations 0] , one usually starts with a known particular solution of a Riccati 
equation without linear term 

f + /,^ = V,W, (6) 

that we call the bosonic Riccati equation. We are not interested here in a so-called factor- 
ization constant that can be placed at the right hand side. The function Vi is an exactly 
solvable potential for the Schrodinger equation at zero energy 

- fpjWb = —^ - Vi{x)wb = . (7) 

The particular solution wi, is usually called a bosonic zero mode. It is connected to the 
Riccati solution through fp = —dwb/dx. Next, changing the sign of the ffist derivative in 
Eq. (H) one calculates the outcome V2{x) using the same Riccati solution 

-f + /,^^V^.W. (8) 

We call the latter equation the fermionic Riccati equation. The function V2{x) is known 
as the supersymmetric partner of the initial potential Vi{x). The corresponding zero mode 
fulfills 





= ^ - V2ix)w, = . (9) 

On the other hand, aheady in 1984, Mielnik [Q studied the ambiguity of the factorization 
of the Schrodinger equation for the oscillator that led him to the general (one-parameter- 
dependent) solution of the companion Riccati equation for that case. In other words, one 
looks for the general solution p| 



dx ■'' dx ■'P 



+ f; = -^ + fi = V2ix). (10) 



The latter equation can be solved for fg by employing the Bernoulli ansatz fg{x) = fp{x) — 
-i^-, where v{x) is an unknown function [^. One obtains for the function v{x) the following 
Bernoulli equation 

^""^"^^ + 2v{x) U{x) = 1 , (11) 



dx 



that has the solution 



T ofc(x) + X 
wl{x 



v{x) = t^^!^^,L^ , (12) 



where 2ob(x) = j^ wliy) dy, and we consider A as a positive integration constant that is 
employed as a free parameter. 

Thus, the general fermionic Riccati solution is a one-parameter function fg{x; A) of the 
following form 

fg{x-X) = U{x) - ±\\n{l,,{x) + \)] = -^[in (-^^^M_\ 1 = ±\^nj,{x-X) (13) 



dx 



dx 



,[Xob(x) + A] 



dx 



where wi,{x]X) = x^{x)+\ - '^^^ range of the A parameter is conditioned by Tq{x) + X ^ 
in order to avoid singularities. This is a well-known restriction ||^. According to the 
super symmetric construction, one can use this general fermionic Riccati solution to calculate 
a one-parameter family of bosonic potentials as follows 



dx ■'' 



+ fa = Vl,, , (14) 



where 



enters the linear equation 



Vi,g = Vi-2-^ In (Tofe(x) + A) (15) 



dbc ~ ^V ' " d^^ Vi^gM^-^ A) = . (16) 

In the limit A — * oo, Eq. (|l^) goes into Eq. (|^) because fg -^ fp and Vi^g — > Vi. 




One can think of Eq. (H) as a generalization of the thermodynamic Riccati equation (^. 
Of interest are the one-parameter oscillator actions fg rather than the 'potentials' Vi^g. Since 
in supersymmetric quantum mechanics Vi^g are the general Darboux-transformed potentials, 
we shall call the fg as Darboux-transformed actions. 




A. The Planck case 

Using / = w' jw in Eq. (|), where the prime denotes the derivative with respect to x, 
leads to the hnear equation 

(17) 

having the particular zero-mode solution Wa = WaSinh(|a;). Thus, the general Riccati 
solution is a one-parameter function fgp{x; A) of the following form 

fgp{x; A) = fp{x) - — [ln(Xoa(x) + A) 




}l0aix)+X], 

Accordingly, the two-variable entropy will also become a parameter-dependent function 

rx 

Sg{x, y; A) = xfgp{y; A) - / fgp{z; \)dz , (19) 



where from all the basic calculations as performed by ACP can be easily repeated. For 
example, to calculate Carnot efficiencies one can use the generalized ACP formula 

1 Tcoid Sg{b, a; A) - Sg{u; A) . 

VC,g = i - 7f 7 TT 7 TT , l^Uj 

Thot Sg{a;X) - Sg{v,u;X) 

and the same values of the parameters as in ACP, i.e., Tcoid = 1/4, Thot = I, a = Phot^i = 1, 
b = Pcoid'^2 = 4:, u = Pcoid^3 = 2c, and v = Phot^A = 2, keeping c as a free parameter. The 
results of this subsection are illustrated in the plots of Figs. (1) - (4). 



B. The vacuum case 

For this case, the particular Riccati solution is the vacuum action fp = fv = |- The 
corresponding zero mode is wy oc e'*^/^. The usual entropy sy(x, x) of the vacuum fluctua- 
tions is zero as a result of a simple calculation, whereas the modifled entropy has a kink-like 
behavior between the | (bosonic) solution and the — | (fermionic) solution. Plots of this 
case are displayed in Figs. (5) - (8). 



C. The symmetric zero mode: Fermi-Dirac action at negative T 

The cases in A and B could be considered particular cases of the general zero-mode 
Wg = Ae^^l"^ -|- Be'^^^"^ that can be also used as solution in Eq. (|^. The Planck action 
corresponds to A = —B = ^ (antisymmetric zero-mode), while the vacuum case to A = 
arbitrary and B = 0. One can use any other type of zero-modes. For example, the symmetric 



zero-mode Ws = WsCosh.{^x) (see Figs. (9) - (10)) is an interesting case since if we trace 
back to the action we get 

^'^~2^ expi-hx) + 1 ' ^^^^ 

i.e., a Fermi-Dirac action for negative x. We have shown in a previous paper that the 
A parameter is equivalent to the quotient A/B 0. Thus, the general Riccati solution 
introduces effects of the second linear independent solution. The problem then turns into 
a subtle interpretation of the mathematical results. We have found at least one possible 
physical significance. For the Planck case, the second linear independent zero-mode is exactly 
the cosh function and the corresponding action is the aforementioned Fermi-Dirac action of 
negative x. We attach now the minus sign to the temperature parameter in x = u/T and 
recall that the issue of negative absolute (spin) temperatures first appeared in Physics in 
1951 when Purcell and Pound were able to produce sudden reversals of the direction of an 
external magnetic field applied to a crystal of LiF 0]. Since then many other experiments 
with negative temperatures have been devised in nuclear spin systems and the corresponding 
'violations' of the second law of thermodynamics were a subject of discussion ^ . Thus, the 
one-parameter Darboux transformations of the Planck action are a way of introducing upon 
it the effect of a Fermi-Dirac action of negative absolute temperature. On the other hand, 
in the figures (9) - (10) we introduce effects of the Planck action of positive temperature on 
the Fermi-Dirac action of negative temperature. 



III. THE DARBOUX GENERALIZATION OF THE 
FLUCTUATION-DISSIPATION THEOREM 

Another interesting application refers to dissipative RLC systems where the oscillator 
action enters Nyquist- Johnson spectral power noise distributions of the type (fluctuation- 
dissipation theorem 0) 

P{u;,P) = -R{u;,(3)fp. (22) 

One can think of the corresponding generalization 

P{u;,(3-X) = -R{u;,(3)fg{(3u;-X) (23) 

vr 

and hope to study even in simple experiments the significance of the present approach 
predicting a Darboux generalization of the fluctuation-dissipation theorem. 



IV. CONCLUSION 

In this work, Planck's thermodynamic oscillator action is generalized to a one parameter 
Darboux family of actions. We also consider the bosonic vacuum case separately in the same 
way. The Planck action and the pure vacuum case correspond to the asymptotic limit of 



the Darboux parameter A — > oo. In the Planck case, all the other A cases correspond to 
a system made of bosons at temperature T interacting with an equal system of fermions 
at temperature — T. In the vacuum case, the A 7^ oo cases describe the interaction of the 
bosonic and fermionic vacua. The efficiencies of the ideal oscillator Carnot cycles based on 
the Darboux-modified Planck and vacuum entropies are calculated along the lines described 
by Arnaud, Chusseau and Philippe. Systems of negative Kelvin temperatures are hotter than 
those of positive T Q and therefore they always represent the hot bath. In real, dissipative 
cases, the same type of generalization is suggested for the fluctuation-dissipation theorem. 
Finally, we mention that a multiple-parameter Darboux generalization is also possible |T^ . 
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Fig. 1. The one-parameter Darboux-modified Planck action fg{x; A) as a function of x and A. 




Fig. 2. The standard and the one-parameter entropy functions. There are only small differences 
between them. For more details see Fig. 3. 




Fig. 3. The difference between entropies As = s„ — sp. 




Fig. 4. Darboux-modified Carnot efficiencies for tlie Planck case. The plane at the height 0.75 
corresponds to the maximum Carnot efficiency for the parameters used by ACP. 




Fig. 5. The Darboux-modified zero-point action. It starts at the normal ^ value at a; = and 
goes to —2 at large values of x. The shape is that of a usual kink (switching) function between 
±2 for large values of the parameter. For small A values see Fig. 6. 




Fig. 6. The same as in the previous figure but for a different orientation to emphasize that at low 
values of A the general Riccati thermodynamic vacuum kink deviates from the common definition 
of a kink and even turns singular for \c= {x — sinha;)/2. 
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Fig. 7. The usual vacuum entropy and the Darboux-modified one. 
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Fig. 8. The same as in the previous figure but for a different orientation. 




Fig. 9. The one-parameter Darboux-modified Fermi-Dirac action of negative T. 




Fig. 10. The one-parameter Darboux-modified Fermi-Dirac entropy of negative T. 



Plots not discussed in the text of the work. 
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Fig. 11. Heat capacity for the one-parameter Planck case compared to the standard case. 




Fig. 12. Heat capacity for the one-parameter vacuum case compared to the standard case. 
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Fig. 13. Heat capacity for the one-parameter Fermi-Dirac case of negative T compared with the 
nonparametric case. 



